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Abstract 

A facilit y 311(1 apparatus are described which determine stiffness, damping, and added-mass rotordynamic 
coefficients plus steady-state operating characteristics of high speed hydrostatic journal bearings. The apparatus has 
a current top speed of 29800 rpm with a bearing diameter of 7.62 cm (3 in.). Purified warm water, 55 °C (130 °F), 
is used as a test fluid to achieve elevated Reynolds numbers during operation. The test-fluid pump yields a bearing 
maximum inlet pressure of 6.9 Mpa (1000 psi). Static load on the bearing is independently controlled and measured 
Orthogonally mounted external shakers are used to excite the test stator in the direction of, and perpendicular to, the 
static load. The apparatus can independently calculate all rotordynamic coefficients at a given operating condition. 

Introduction 

Experience with the SSME (Space Shuttle Main Engine) has demonstrated definite limits for the successful 
operation of ball bearings in liquid oxygen (LOX). Currently, bearings in the HPOTP (High Pressure Oxygen 
Turbopump) of the SSME experience accelerated wear at full power level. The balls in these bearings simply get 
smaller, rapidly. The problem has proven largely intractable. As a result, hydrostatic bearings have been proposed 
lor many new turbopump applications because of their long lifetime, low friction factors, low wear, and their ability 
to use low-viscosity lubricants. For cryogenic applications, the bearings will be pressurized from pump-discharge 
now. At zero speeds the bearings will be flooded but unpressurized. Rubbing will be experienced at start up and 
shut down as investigated by Scharrer et al. (1992 a, 1992 b). The present research concerns steady-state fully 
pressurized bearing operation. 

The test apparatus and facility described here is used to develop experimentally validated tools to predict 
steady-state hydrostatic bearing operational data (resistance torque, static load, flowrate, temperature and pressure 
distributions, etc.) and rotordynamic coefficients for vibration analysis. By rotordynamic coefficients, we refer to the 

stiffness K, damping C, and added-mass M coefficients which are used in the following linearized force-displacement 
model for bearings 



Ajc ^ C x * C xy fAil M xy Uj? 

41 K» W + k cj Wj + [M yx M yy \ (Ajf 


( 1 ) 


Aj,)define ^ motion of * e bearing rotor relative to its stator, and are the components of the 

fluid film reaction force acting on the rotor. Identification of the rotordynamic coefficients of Eq.(l) is a central 
objective of this research project 

This paper describes the test-facility design requirements necessary to identify the rotordynamic coefficients 
and presents a sample of the resulting data for one bearing. The contents include results for a 7.62 cm (3 in.) 
diameter, square-recess, smooth-land, annular-fed and orifice-compensated hydrostatic journal bearing with an L/D 
rauo of 1 and a C/R ratio of 0.003. The design requirements and analysis techniques for determining the steady-state 
operational data are detailed thoroughly by Kurtin, et al.(1991) and will not be repeated here. Subsequent publications 
will provide test results for hydrostatic bearings with assorted recess shapes, distinct land patterns and roughness 
varied area ratios, and tangential-injection against rotation. 


The research reported here was supported in part by Rocketdyne, Division of Rockwell International 
under Contract No. R94QBA89-032182, NASA through the Texas A&M University Center for Space 
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Apparatus and Facility - An Introduction and Overview 

The Hydrostatic Bearing Test Stand illustrated in Figure 1 consists of the following major components. A 
base, fabricated from welded steel plates, supports the test section, gearbox, and drive motor. In the test section, two 

stainless steel pedestals spaced about 38 cm (15 
inches) apart support a 7.62 cm (3 in.) diameter 
stainless steel high-speed rotor. The rotor is driven 
by a 93.2 Kw (125 HP) variable-speed electric 
motor through a 7:1 ratio speed-increasing gearbox 
via a disc-type coupling. The test-bearing stator 
supports itself on the rotor at a position midway 
between the support pedestals. A pneumatic 
cylinder provides a fixed-direction side load to the 
bearing. 

Two orthogonally mounted hydraulic 
shakers are attached to the stator through isolation 
stingers as illustrated in Figure 2. The shakers 
provide excitation parallel and perpendicular to the 
static side load along the X and Y directions, 
respectively. Each shaker can excite the stator with 
measured dynamic loads to 4450 N (1000 Lb) at 
frequencies to 1000 hz. The stingers were designed 
following guidelines by Mitchell and Elliott 
(1984). The static and dynamic loading 
arrangement of Figure 1 and Figure 2 is most 
closely related to that of Kanki and Kawakami 
(1984) who use orthogonal static loads from 
pneumatic cylinders in series with electromechanical shakers to excite the stator in a liquid-seal test rig. Glienicke 
(1966-67) is normally credited with the initial development of this rig type. Tonnesen and Hansen (1982) used this 
type of hydrostatic support bearings setup in static tests to investigate temperature effects of hydrodynamic bearings. 

Preloaded axial cables arranged in an opposing 
trihedral configuration (Figure 3) constrain the shaft during 
stfliir loading and dynamic excitation. A pitch mode 
instability of the stator with rotor onset speed around 17,000 
rpm was also eliminated by these cables. This pitch restraint 
concept came from Nordmann and Massman (1984), who 
used a similar design on a liquid-seal test rig. The pitch- 
restraint cables have an additional advantage in insuring that 
the stator remains parallel to the rotor during excitation. 

Alternatively, they have been used to deliberately misalign 
a test bearing during excitation. 

Reluming to Figure 2, note that accelerometers are 
provided to measure the stator acceleration in the X and Y 
directions. Although not illustrated in Figure 2, the relative 
deflections Ax(t),Ay(t) between the bearing rotor and stator 
are also measured with high-resolution eddy-current 
transducers. In summary, dynamic measurements at the 
bearing for rotordynamic coefficient identification include 
the dynamic excitation forces, the stator acceleration, and 
the relative rotor-stator deflections in the X and Y 
directions. 

Purified water at 55°C (130°F) is supplied to the 



Figure 2 End view of a test bearing 
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support and test bearings via a centrifugal pump with pressure and flow rate capacities of 6895 Kpa (1100 psi) and 
303 liters/min (80 gpm), respectively. The main pump supplies the water to the lest section through 10-micron 
absolute filters. The test-bearing flowrate is monitored by a turbine flowmeter and adjusted by a remotely-controlled 
valve. Nitrogen-backed accumulators provide a pressurized bearing coastdown in the event of a main pump failure. 
The water is elevated to test temperature by heat input from 


the main pump and test section. Two PID-controlled valves 
mix cooler water from a large 15,000 liter (4000 gallon) 
reservoir to maintain a constant fluid temperature at the test 
bearing inlet. The equipment is protected by an extensive 
logic-controlled monitoring and fault detection network 
capable of partially or totally unassisted safe shutdowns. 

The test apparatus has been designed to facilitate 
rapid testing of multiple bearing configurations. The 
coupling design of Figure 1 allows rapid removal of a test 
rotor, followed by speedy replacement of a lest bearing. A 
complete bearing test is accomplished in one week, 
including bearing and rotor installation, calibration, testing, 
and a computer-generated test report. Actual testing is 
computer driven and takes on the order of one-half day. 



Rotordynamic-Coefficient Identification 

The preceding section introduced the test apparatus and facility and explained the dynamic measurements 
taken at the bearing. The present section explains how the measurements are used to identify the rotordynamic 
coefficients. 

The equations of motion for the stator mass M, of Figure 2 can be written 


M. 



( 2 ) 


Where X„y, are the (measured) components of the stator’s acceleration,//,/, are the (measured) components of the 
input excitation forces, and f Mi , f yh are the bearing-reaction force components. The x and y subscripts in these 
equations identify the X and Y directions of Figure 2. Substituting from Eq.(l) and rearranging Eq.(2) gives 
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Since the mass M, is known, the left-hand vector of Eq.(3) is a known function of time. On the right-hand side, Ax(t) 
and Ay(t) are measured functions of time. The rotordynamic coefficients are determined in the frequency domain 
via the Fourier Transform ST which yields, 
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where 


FfWfo A k =jr(fy, D k =&{ A*) 


(5) 


The elements of the frequency-response function (FRF) H are related to the coefficients defined in Eq. (3) by, 

^xy ~ (K X y~V> 2 M xy ) +/(«C ly ) 

H yx = (K yx - »'M yx ) + y(«c„) 

H yy « (Ay y - a 2 M yy ) + ;'(«C y> ) 

where 0) is the excitation frequency and / = y-1 . 

However, Eq.(4) only provides two complex equations for the four unknown HJ s. To get four independent 
equations, two independent orthogonal dynamic loads are applied to the stator by alternately exciting along the X 
and Y directions of Figure 2. Each of these "pseudo-random" dynamic loads contains 41 sinusoids. The superposition 
of these sinusoids is optimized to provide a composite loading that has a high ^ spectral-line energy to crest-factor 
ratio. The loadings are periodic and have energy only at 10 Hz frequency increments in the range from 40 to 440 
Hz. The two independent loadings result in the four complex equations, 


F xx - M Axx F xy " M Axy 


H xx H xy 

\D D 1 

xx xy 

F yx M f A yx F yy MjA yy 


H „ H yy 

D D 

y* yy] 


where the subscripts ij of the force and displacement matrices correspond to response in the / direction due to 
dynamic loading in the j direction (Figure 2). Appendix A provides a complete discussion of the power-spectral 
density method used to calculate the H u functions. 

Steady >State Operating Conditions 

The test equipment varies the test bearing steady-state operating conditions by independently controlling the 
following four parameters: 

1) Rotor speed, 

2) Test fluid supply pressure, 

3) Static bearing load, and 

4) Test fluid supply temperature. 

The normal operating-condition test matrix includes three running speeds from 10,000 to 25,000 rpm, three supply 
pressures from 4 to 7 Mpa (600 to 1000 psi), and six eccentricity ratios (determined by the static load) ranging from 
concentric to 0.5. Tests are conducted at a single temperature, usually 55°C (130 °F); however, the test apparatus 
has been operated to 29,800 rpm at 71°C (160°F). 

Reynolds-numbers estimates for the test seal at 28,000 rpm, C,/R=.0033, with 160°F water is 


1* = * w C,p/p = 37,400, 
R = /r5+R 2 ‘ * 42,700, 


JL z = 2WC r pl\i = 20,600 

iyR t = 1.82 


( 8 ) 
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( 9 ) 


For comparison, a typical liquid-oxygen lurbopump bearing would have 

*e = 91.300, R, = 113,200, R = 145,400, R e /R j = .81, 

and a typical hydrogen bearing would give 

= 121.300, H t = 366,600, R = 386,200, R^ = . 33 . 


( 10 ) 


R for the LOX and H 2 bearings 
is seen to be about 3.5 and 9 
times higher, respectively, than 
corresponding water values. 

Test Bearing 

As reviewed earlier, the 
test apparatus is used to identify 
rotordynamic coefficients of 
hydrostatic bearings with a 
variety of geometric 
configurations. The results 
presented here are for the 5- 
recess bearing of Figure 4, 
whose physical properties and 
dimensions are described in table 
1. The bearing area ratio is the 
recess area divided by the total 
bearing area. The bearing’s outer 
and inner shells are made, 
respectively, from aluminum and 
bronze. Shallow recess depths 
are typically required in 
cryogenic bearings to avoid 
"water hammer" by minimizing 
the recess volume. Because of 
the small depths, and the 
requirement of a smooth recess 
surface, recesses were 



manufactured by an Electrical 
Discharge Machining (EDM) process. 


. „ nl h Waler enters u 38 llIustrated at ^ l °P l e ft-hand element of figure 4, circulates through a 360* annulus 
and enters the recesses through orifices. Close-coupled, strain-gauge pressure transducers 

dynamic pressures within each bearing recess. The unwrapped faring view at teSmtflZt 4 Ld Te end 

respectively^r^ measurements made on a test bearing, with ©, 0, 0, 0 and ® denting. 

' PreSSUreS ’ “ ,P “ atUrcS - and “ etali0nS ' 2 ^ 


Figure 5 provides a detailed view 
access and easy replacement of each orifice, 
from the orifice exit to the recess. 


of the orifice and recess cross section. The apparatus design provides 
bserve that there is a 90 expansion flow cone for a transition section 
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Table 1. Test-bearing physical parameters (Refer to figure 4 for dimension definitions.) 


Mass 

Diameter (O.D.) 

Diameter (I.D.) 

Radial Clearance 
Length (L) 

No. of Recesses 

Recess dimension A 

Recess dimension B 

Recess Volume 

Area Ratio 

Recess depth 

Orifice Diameter 

Orifice Supply Line Volume 

Land Roughness (peak-peak) 


11.34 kg (25 lb) 

76.441 mm (3.0095 in) 

76.200 mm (3.0000 in) 

0.12 mm (0.0048 in) 

76.2 mm (3 in) 

5 

27.0 mm (1.064 in) 

27.0 mm (1.064 in) 

0.184996E-6 m 3 (0.0112891 in 3 ) 

0.2 

0.259 mm (0.010 in) 

2.49 mm (0.098 in) 

0.128994E-6 m 3 (0.00787173 in 3 ) 
0.330 pm (13 pin) 


Table 2. Mnemonic definitions for measurements in figure 4. 


measurement 

MNEMONIC 

LOCATION/DIRECTION 

RELATIVE MOTION 

Ml 

FRONT / X- 


M2 

FRONT / Y- 


M3 

REAR / X 


M4 

REAR/ Y 


M5 

REAR / X- 


M6 

REAR / Y- 

PRESSURE 

PI 

RECESS 1 


P2 

RECESS 2 


P3 

RECESS 3 


P4 

RECESS 4 

— 

P5 

RECESS 5 

" 

P6 

INLET ANNULUS 


P7 

BACKSIDE ANNULUS 

TEMPERATURE 

T1 

ANNULUS 


T2 

LAND AT RECESSES 


T3 

RECESS 1 


T4 

LAND AT EXIT 

acceleration 

Ax 

CENTER /X 


Ay 

CENTER / Y 

FORCE 

Fx 

CENTER / X 


Fy 

CENTER / Y 
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Representative Dynamic Test Data 
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Recall from Eq.(6) that K„ and M„ coine from the real pari of H„, and C„ is obtained from 
the imaginary part of H xx , etc. Figure 7 illustrates the rotordynamic coefficients versus eccentricity ratio 
for a running speed of 24,600 rpm, a supply temperature of 54.4'C (130’F), and a range of supply 
pressures. The rotordynamic coefficients are seen to be generally insensitive to eccentricity ratio. Recall 
from figure 4 that the bearing is loaded toward the center of 
a recess, and the bearing is not symmetric. Hence, the direct 
coefficients (e.g. and M^) are not necessarily equal at 
zero eccentricity. Note the large value for direct and cross- 
coupled mass coefficients which arise because of the elevated 
Reynolds number. 

Figure 8 illustrates whirl-frequency-ratio f„ versus 
eccentricity ratio for three supply pressures at 24,600 rpm. 

The whirl-frequency ratio definition developed by San 
Andres (1991b) was used in these plots and depends on all 
the rotordynamic coefficients. Although not illustrated, f a 
increases slightly with decreasing running speed. In the 
absence of other rotordynamic forces, the onset speed of 
instability of a flexible rotor supported by plain journal 
bearings is approximately equal to the first critical speed 
divided by the whirl-frequency ratio, Lund (1966). 

Unfortunately, these test results show that f w s 0.5 for the 
hydrostatic bearings which is no better than a plain journal 
bearing. Constaninescu (1969) predicts this outcome for Figure 8 Whirl-frequency ratio 
hybrid gas bearings. versus eccentricity ratio for three 

supply pressures 

Uncertainty Analysis 

A "dry-shake" test (one in which the bearing fluid film consists only of air at atmospheric 
pressure) was conducted to examine the magnitude of the stiffness and damping imposed on the stator 
by the cable, stinger, and other attachments. These values are negligible compared to the stiffness and 
damping of all bearings tested so far. 

Table 3 contains estimates of the upper limit bias errors the precision indices and the 
absolute uncertainties U RSS of the coefficients for the test case of Figure 6. These are obtained by 
calculating bias and precision errors for each element of the force and displacement matrices of Eq.(7) 
then numerically propagating these errors through the functional relationship between the elements and 
the resulting coefficients. The method used conforms to the ANSI/ASME standard on Measurement 
Uncertainty (1986). The coefficients fall within VALUE +/- U RSS for 95% coverage. 



Table 3. Measurement Uncertainties for co = 24600 rpm, zero eccentricity, and 6.9 MPa (1000 psi). 


COEFFICIENT 

VALUE 

Brss 

Srgs 

Urss 

UNITS 

K xx 

146 

2.82 

1.82 

4.60 

MN/m 

K X y 

136 

2.21 

1.34 

3.47 

MN/m 

Kyx 

-133 

2.11 

1.86 

4.27 

MN/m 

Kyy 

132 

2.46 

1.30 

3.59 

MN/m 

C X x 

112 

2.26 

0.74 

2.71 

kN/m 

C X Y 

14.4 

2.60 

1.42 

3.85 

kN/m 

Cyx 

-11.3 

2.39 

0.69 

2.76 

kN/m 

Cyy 

105 

2.16 

1.54 

3.77 

kN/m 

M XX 

8.8 

1.06 

0.99 

2.24 

kg 

M xy 

-5.1 

0.50 

0.79 

1.66 

kg 

Myx 

2.2 

0.44 

1.00 

2.04 

kg 

My Y 

2.9 

1.06 

0.59 

1.58 

kg 
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Summary, Conclusions, and Comments 

The results and discussions included here present and explain the design of a new test system 
for hybrid (combined hydrostatic and hydrodynamic) bearings. The tester can measure static parameters 
including pressures, flowrate, temperatures, drag torque, etc plus a full complement of force/displacement 
rotordynamic coefficients. The test system expedites rapid and flexible testing of fluid-film bearings. 
It combines proven design elements of earlier experimenters and state-of-the-art excitation and 
identification techniques to deliver high-quality data. 

Space limitations have forced us to present only a very limited selection of the data which have 
been developed for a single "generic" bearing. Tests have been completed for approximately 30 bearings 
as of 1 August 1992. 

As noted in the introduction, this test program aims to produce test data to help in developing 
and validating predictions which have been under concurrent development by San Andres (1990a, 
1990b). San Andres uses a 2-D, bulk-flow Navier-Stokes model for flow in the bearing land, a 
compressible flow continuity requirement for the recess, a discharge-coefficient model for the orifice, 
and concentrated pressure drops for flow leaving the recess and entering the lands. The analysis applies 
for rectangular recesses and radial flow injection. Comparisons between measurements and predictions 
are generally good for rectangular-recess bearings and are presently being prepared for publication. 

A central question involved in the development of the present apparatus concerned the required 
excitation device. Specifically, should impulse devices or shakers be employed? Impulse excitation is 
a standard approach for modal identification, Ewins (1986), and was pioneered by Nordmann for bearing 
parameter identification and subsequently for annular seals, Nordmann and Massman (1984). Hydraulic 
shakers were first used for parameter identification of annular seals by lino and Kaneko (1980). 

Price is a central advantage of impulse devices over hydraulic shakers. The impulse "guns” 
which were initially purchased for this project cost less than one tenth the eventual cost of the hydraulic 
shakers. The decision to shift from impulse to shaker excitation, despite financial considerations, 
involved the quality and volume of data to be acquired. The hybrid bearings were extremely noisy. 
Pressure pulsations in their recesses generated unwanted relative bearing-rotor amplitudes of the same 
magnitude as that created by external excitation. As detailed by Rouvas et al. (1992), impulse excitation 
can be made to work, but only with many repetitions and averaging in the time domain. However, the 
FRF data developed by Rouvas et al. with impulse excitation never approaches the high quality shown 
in figure 6, despite averaging over 200 separate impulse excitations. 

An additional practical problem with impulse excitation concerned tip selection for the impulse 
guns. Soft tips were required to deliver the required low-frequency force excitation, but they would not 
have survived testing a bearing over a large test matrix with 200-plus excitations per data point. A 
parallel identification concern involved the comparatively large peak rotor/stator amplitudes developed 
by impulse excitation which were necessary to give adequate amplitudes at discrete frequencies. 

Having acquired and installed the hydraulic shakers, the nature of the excitation and 
identification procedures remained to be resolved. Pseudo-random excitation is used because it is fast, 
has no leakage in the analysis. The spectrum is shaped to only excite frequencies in the range of interest 
with energy only at each DFT analyzed spectral line, and only a few averages are needed due to its low 
crest factor and good signal-lo-noise ratio (Herlufsen, 1985). This excitation technique, along with 
identification procedures as outlined by Bcndal and Piersol (1986) and extended by Rouvas and Childs 
(1992)“ eventually achieved the high-quality FRF’s of figure 6.. 

The lest results developed here show one principal limitation to the planned use of hybrid 
bearings in cryogenic turbopumps, specifically, a whirl-frequency ratio near 0.5. With this whirl- 
frequency ratio, running speeds above twice the rotor’s first critical speed are likely to yield unstable 
motion. Hence, recent developments within this program have considered the following options for 
improving rotordynamic stability by reducing f„: 


4 See Appendix A 
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il!^ eC ! ed injecti0n agalnst rotation - Tondl ( 19 67) initially proposed injecting the fluid 
against rotor rotation as compared to the radial injection of figure 4. 

(,982) W— —» for a^iar 

“PP|‘ oach ® s acl 10 reduce the average tangential velocity within the bearing. Favorable 
C ° mPlCle<l b °' h ‘ yPCS ° f impro, ' d bc " ri "« ronfiguracion* and will be 

turhoniimm e H mPhaS1S ? ^ W ° rk ^ bcen on lhe application of hybrid bearings to cryogenic 
nnS 7 ■ e ^’ ^ ey presem 0bvi0US alternatives as product-lubricated bearings for commS 

Sv^r and0kada < 1992 > hav * recend y reported the results of a 6-year development p7mm 
Tokyo Gas Co. implementing hybrid bearings into vertical liquified-nalural-gas pumps. & 

• t , A no ^ by KurUn et al - ( 1991 ). two practical problems have been encountered which might 
restrict the applications of hybrid bearings to commercial pumps. First, considerable cavitation erosfon 

H V ™ d part 18 made from 174PH stainless steel with a high yieldTtrength 
Tes^ of a modified bearing with the orifice exit flush to the recess* showed no caviLTdaS on 

h Ut S 7 e r ° Ughening of the rotor - Given that the recess pressure is on the order 
of 3.5 MPa (500 psi) this damage could be resulting from air bubbles coming out of solution in the 
water versus an actud phase change of the water. Long-term testing to investigate ^7^11 
damage is recommended before these bearings are selected for commercial pumps. Tests with the product 
fluid which will actually be used are advisable, since almost any fluid will be kinder from a cavitation 
viewpoint than water. Rocketdyne personnel report zero damage with liquid hylgen 

nriainaii ^ 7°^ Pr ° b Cm C0nCQms dama « e t0 the faring due to fluid-borne debris. Our flow loop 
strim^ V^i 3 °' micron " ominaJ fi,ter which wa s inadequate in removing particles from the flow 
su-eam Particles entering the bearing would lodge at the comers of the recesses and then create 
circumferential damage grooves around the bronze stators. Switching to a 10-micron absolute filter 

thi^Sg ty^ N ° ne ^ 1CSS * CarC iS suggested in providing ade( l uale ^stream filtration for 


5 T he divcrging transition section from the orifice exit to the recess is eliminated. 
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Appendix A. Spectral-Density-Averaging for Rotordynamic-Coefficient 
Identification 

Spectral density averaging technitjues are used to obtain the frequency-response functions (FRF) 
of Eq.(7), which defines a multi-input multi-output (M1MO) system. The inputs are D xx , D xr D )X , and 
D and the outputs are F„ - F Xf - M,A xr F yx - Af/4^, and F- M/i„. 

Before calculating the FRF for this MIMO system, first consider the best solution for the FRF 
of a single-inpul single-output (SISO) system. A S1SO system with a clean input a and a noise 
contaminated output b will be examined. According to Bendal and Peirsol, (1986) an unbiased estimate 
H, of the FRF for this type of SISO system is given by. 


H - 



(A.I) 


with G » being the cross spectral density and G„ the auto spectral density. For discrete data, (which we 
have) an FFT algorithm is used to calculate the Discrete Fourier Transforms (DFT’s) of Eq.(A.l). In 
particular, the spectral densities are calculated by, 




(A.2) 


where n d is the number of statistically independent blocks of data, N is the number of data points in each 
data block. A/ is the time increment between sampling of each data point, A (to) = the complex conjugate 
of ria(t)} and P( to) = 7[p(t)) >s used to represent the DFT here) 

We wish to apply the unbiased estimator for the SISO system described above to the MIMO 
system of Eq. (7). First, recall that we have 32 statistically independent blocks of data (n d = 32), each 
block containing 1024 (N = 1024) data points for each of the dynamic loads, accelerations, and relative 
motions whose J appear in Eq.(7). Next, we identify each independent dynamic load of Eq.(7) (F„ and 
F ) as a clean input a, and each of the resulting elements of the force and displacement matrices as 
noise contaminated outputs b. Eq.(A.l) is then applied separately to each individual element of the force 
and displacement matrices, y ielding an unbiased estimate for each input and output element contained 
in these matrices. The resulting equation, 
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yields improved estimates for the FRF when solved. . . 

Once Eq.(A.3) is solved, the rotordynamic coefficients are extracted from the FRF definitions 
of Eqs.(6) by performing a least-squares curve fit on the real and imaginary parts of each element of the 
FRF matrix. Plots of typical real and imaginary parts for Hxx are shown in Fig. 6. The stiffness, 
damping, and added-mass coefficients of Eqs. (1) and (3) have now been identified. A more detailed 
discussion is provided by Rouvas and Childs (1992). 
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